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Abstract 

We study a zero-range process where the jump rates do not only de- 
pend on the local particle configuration, but also on the size of the system. 
Rigorous results on the equivalence of ensembles are presented, character- 
izing the occurrence of a condensation transition. In contrast to previous 
results, the phase transition is discontinuous and the system exhibits ergod- 
icity breaking and metastable phases. This leads to a richer phase diagram, 
including nonequivalence of ensembles in certain phase regions. The paper 
is motivated by results from granular clustering, where these features have 
been observed experimentally. 

keywords, zero range process; discontinuous phase transition; equivalence of 
ensembles; metastability; ergodicity breaking; granular clustering 

1 Introduction 

The zero-range processes is an interacting particle system introduced in [|29l . 
which has recently attracted attention due to the possibility of a condensation 
transition. A prototype model with space homogeneous jump rates that exhibits 
condensation has been introduced in flU. When the particle density p in the sys- 
tem exceeds a critical value pc, the system phase separates in the thermodynamic 
limit into a homogeneous background with density pc and a condensate, that con- 
tains all the excess particles. This phase transition is by now well understood on 
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a mathematically rigorous level for general zero-range processes [fTTl . and has 
been applied to model clustering phenomena in various fields (see ifTOl and ref- 
erences therein). In one dimension, a mapping to exclusion models gives rise 
to a criterion for non-equilibrium phase separation [19J. Further rigorous results 
on the zero-range process include a proof of condensation even on finite lattices 
[fT2ll . and a refinement of the results in [[TtII . which implies a limit theorem for 
typical density profiles in case of condensation [|24l . Regarding the background 
density as order parameter, it has been shown in a general context (including dif- 
ferent particle species) that in spatially homogeneous zero-range processes with 
a stationary product measure condensation is always a continuous phase transi- 
tion [fT6l . Recently, investigations have been further extended to open boundaries 
where particles are injected and extracted [21] and heuristically to various gener- 
alized models. Those include a non-conserving zero-range process that exhibits 
generic critical phases f2\, zero-range processes with non-monotonic jump rates 
leading to multiple condensate sites [28 J, or mass transport models with pair- 
factorised stationary measures that give rise to a spatially extended condensate 

m. 

In this paper we study the condensation transition in a generalized zero-range 
process where the jump rates depend on the system size. The motivation for 
this study comes from experiments on granular media reported in [|26l [34l l32l . 
Granular particles are distributed uniformly in a container which is divided in 
several compartments. When shaking the container, the particles start clustering 
in some of the compartments and after equilibration, almost all particles form 
a "condensate" in one of the compartments. The phenomenon is robust for a 
variety of shaking strengths and a gas-kinetic approach lead to a simplified model 
equivalent to a zero-range process where the hopping rates depend on the number 
of compartments [|71 [34l [32l [33l. In an alternative activated-process approach it 
can be modeled by a zero-range type process, where the jump rates depend on 
the total number of particles in the system [2T,'4'|, and both approaches have been 
summarized in [30|. A heuristic analysis of the behaviour of the order parameter 
agrees with experimental observations and shows that generically the transition is 
discontinuous and the system exhibits hysteresis and metastability. This analysis 
suggests that the discontinuity is due to the dependence of the jump rates on the 
total number of particles or the number of compartments, respectively. 

To treat this phase transition on a rigorous level, we present a detailed analysis 
of a simple prototype model with system-size dependent jump rates, for which we 
derive results in the context of the equivalence of ensembles analogous to [[T7l[T6l . 
From a mathematical viewpoint our system provides an interesting example, since 
the origin of the phase transition is due to a non-standard behaviour of the grand- 
canonical measures, in particular the lack of a law of large numbers. This leads to 
a richer behaviour than in previous models, which can be fully understood only 
by studying the canonical measures as well, which is not the case for zero-range 
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processes with fixed jump rates [fT6ll . The mathematical structure is also different 
from standard results on systems with bounded Hamiltonians JH |3T1 . We also 
show how our findings can be directly generalized to a process where the jump 
rates depend on the total number of particles, rather than the size of the lattice. 
To establish the link between the stationary distribution and dynamics we include 
a discussion of metastability and the life times of metastable phases, which are 
compared to Monte Carlo simulation data. Our results can be generalized heuris- 
tically to a large class of systems, including models of granular clustering, as is 
explained in a forthcoming publication LISJ . 

The paper is organized as follows. In the next section we introduce the model 
and show its phase diagram, which summarizes our results. In Section 3 we study 
canonical and grand-canonical stationary measures and the equivalence of ensem- 
bles is discussed in Section 4. In Section 5 we present results on metastability and 
in Section 6 on the extension to a dependence on the number of particles in the 
system. In the discussion in Section 7 we give a detailed comparison with previ- 
ous results. 



2 Model and results 

We consider a zero-range process on a translation invariant lattice Al of size 
\Al\ = L. The state space is given by the set of all particle configurations, 

XL = {r] = iV'.)xeA, ■■ V. e N} , (1) 

where the number of particles per site can be any non-negative integer number. 
With rate gR{r]x) one particle leaves site x G Al, and jumps to another site y with 
probability p{y — x). To avoid degeneracies, we require the jump probabilities 
{p{x) I X G A^,} to be irreducible and of finite range, i.e. p(x) = if |x| > C 
for some C > 1. Under these conditions our main results are independent of the 
actual choice of p. Since they cover the basic novelties of the paper, we restrict 
ourselves to the jump rates of the form 

9Rik) = l^2 [ktn for^>l' ^?(0)=0, (2) 

where Cq > Ci > 0. The rates are piecewise constant and the location of the jump 
is given by the parameter R > 0, which depends on the system size L, such that 

R ^ oo and R/L ^ a as L ^ oo , (3) 

where a > is a system parameter. The most interesting case we will consider is 
a > 0, but we will also discuss a = which depends on the asymptotic behaviour 
of i? as L tends to oo. The same model has already been mentioned in [9| for 
fixed R. There is no phase transition in this case, but for large R one observes a 
large crossover, i.e. convergence in the thermodynamic limit is very slow. 
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Figure 1: Stationary phase diagram for generic values of cq > ci. The four phases F{E) 
ip < Pc), F {pc < p < Pc + a), F/C {pc + a < p < ptrans) and C/F {p > ptrans) are 
explained in the text. 

Left: Phase diagram in terms of a Q and the particle density p. Right: Background 
density as a function of p for a = 0.5. Full lines are stable, broken lines metastable. 



The generator of the process is given by 

It is defined for all continuous cylinder functions f E C (Xl). Since we define the 
process only on finite lattices, there are no further restrictions on initial conditions 
or the domain of the generator as opposed to zero-range processes on infinite 
lattices (cf. yj). We do not specify the geometry or the dimension of the lattice, 
since our main results on the stationary distribution do not depend on these details. 
The only requirement is that the lattice is translation invariant, or more generally, 
(f)x = const, is the only positive solution to the difference equation 

<l>x= (j)yp{x - y) . (5) 

Note that no particles are created or annihilated and the number of particles is a 
conserved quantity. Under our assumptions on p and g there are no other conser- 
vation laws that would lead to degeneracies in the time evolution. 

For fixed L, also i? is a fixed parameter and known results on stationary mea- 
sures for zero-range processes apply (see e.g. [lOJ and references therein). The 
stationary weight w^{r)) for this process is of product form, 

^kiv) = n 

where the single-site marginal is given by 



=0 
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Here the empty product (for A; = 0) is understood to be unity. 

The results we derive in the following sections are summarized in the station- 
ary phase diagram in Figure [T] in terms of the conserved particle density p and 
the parameter a In the fluid phases F{E), F and F/C the stationary mea- 
sure concentrates on homogeneous configurations with bulk density p. In phase 
F{E) for p < Pc the canonical and grand-canonical ensembles are equivalent (see 
Section 4), and in phase F / C for pc + a < p < ptrans there exists an additional 
metastable condensed phase, which has a lifetime exponential in the system size 
(see Section 5). Typical condensed configurations have a p-independent homoge- 
neous bulk distribution with density pc < p, where the excess particles condense 
on a single lattice site. In phase C/F, i.e. for p > ptrans, the condensed phase 
becomes stable and the corresponding fluid phase metastable. On top of metasta- 
bility the order parameters are discontinuous as a function of the density p, and 
therefore the condensation transition is discontinuous. 



3 Stationary measures 
3.1 Grand-canonical measures 

Since the state space is discrete we will identify measures /u({t7}) with their 
mass functions p(r}) in the following to simplify notation. For each R and L there 
exists a family of stationary product measures z/^^ with single site marginal 

iyiRik) = J^^WRik)<P'. (8) 

The marginal is well defined for fugacities G [0, Ci), since the tail behaviour 
of the stationary weight ^ is Wnik) ~ for all fixed R. The single site 
normalization is given by the partition function 

.HW = f = (9) 

and the expected particle density under the measure v^j^ is given by 

Pr{4>) = {Vx)^i = 09^(log^i?,(0)) = 

.0.fl+i i^+l + 0/(c,-0) 

Note that is strictly increasing in and that for every fixed R, Pr{(P) 00 

as Ci. So for all densities p > there exists 0j?(p) such that the measure 
^^r{p) r density p, i.e. product measures exist for all densities. But the single 
site marginals of these measures still depend on R and therefore on the system 



5 



size L. Since R oo as L oo, the marginal ([8]) converges pointwise to a 
simple geometric distribution, i.e. for each A; G N, 

^iR{k)^^U{k) = ^M/cof with zU4>) = ^^. (11) 

This convergence holds for each fixed < ci, but it is not uniform in 0. The 
limiting product measure u^^^o is defined for all < cq. We denote the particle 
density with respect to this measure by 

Poo(0) := (?7x) 1 = 0<9</,(log2;oo(0)) = ^ , , (12) 

<^,oo Co — (p 

and its inverse is given by 

0oo(p)=Co— (13) 

1 + p 



Since convergence (11) only holds for < ci we define the critical density 

Pc ■= Poo(Cl) = — - — < OO . (14) 
Co — Ci 

Note that with this definition 0oo(Pc) = ci- In the following we summarize some 
straightforward consequences of these definitions. 

Proposition 1 For all (p < ci, v^ j^ i^^^oo weakly or, equivalently, 

^f)'^in^^f^^^^- ^sL^^forallf eCoMl\\^ (15) 
and Pr{4>) Poo(0)- For all p > we have 

Mp) , P<Pe . r , p<Pe 

Cl , P>Pc '^^(^^■^ 1 Z^ci,oo , P>Pc 



where (pn is the inverse of ( [70| ) anJ second convergence holds in the weak 
sense as in (1751). 



Proof. (11) implies pointwise convergence of arbitrary n-point marginals z/^ ^ 
and in general this is equivalent to convergence of expected values of cylinder test 
functions, as long as they are bounded. This does not directly imply convergence 
of the unbounded test function rfx which yields the density, but pR{(j)) Poo(0) 
follows by direct computation from ( [TO] ). 

Since pR{(j)) and its inverse are continuous for < ci or equivalently p < pc, 
we have (f)R{p) — * 0oo(p)- Since 0oo(p) < ci and 2:00 is a continuous function, 
inserting 0i?(p) in ^ yields as L — > 00 

zr{Mp)) = ^oo{Mp))(i+(—Y^' """'7. ') ^^oo(0oo(p)()17) 

V ^ Co / Ci—<Pr{P) J 



' Co,6 (^) denotes the set of all bounded, continuous cylinder functions f : X ^ M.. A cylinder 
function depends only on the particle configuration on a fixed finite number of lattice sites. 
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Therefore, we have pointwise convergence of the marginals as in (jlll) and '^^^(p) ji - 
^(t>oo{p),oo weakly for p < Pc analogous to above. For p > pcto leading order 

0i?(p)-ci- — , , , , ^= ci asL^cx). (18) 

For p = Pc the correction has a different power (^)^^^ which leads to the same 
behaviour as for p > pc. Inserting in (|9]) this yields analogous to ( 17 ) 

MMp)) = ^oo(0«(p)) (^1 + i^Y^'^j^) ^ ^-(ci) L 

so that Jy^^^p)^ji Jyci,oc weakly. □ 
Note that by Proposition [1] the density does not converge if p > pc since 

Pi?(0i?(p)) = P 7^ Pc = Poo(ci) , (20) 

and the variance of r]x even diverges as 

/Cg\iJ/2-l 

Var{r]^) = (f)d^pR{(j))\^^^^^^^ - '^^^'^J ^^^^ 



Therefore there is no standard law of large numbers for the measures 
when p > pc- In particular one can show the following. 



<t>R{p),R 



Proposition 2 For each L let r]^ , . . . ,r]j^ be iid random variables with distribu- 
tion '^^^(^p) ji and assume that R ^ log L. Then 

tV^t^I^ .P<Pc almost surely . (22) 
tl^ I Pc , P > Pc 

Proof, see appendix 

Note that for p > pc(2Q) holds due to very large values r]^ ~ (f^)^^^ having very 

/? /9 I 

small probabilities (^) , which also leads to divergence of the variance (21 1. In 
turn, the small probabilities lead to almost sure convergence of the sample mean 
to Pc < p, which is a non-standard strong law of large numbers. The breakdown of 
the standard strong law coincides with the region of nonequivalence of ensembles, 
as has been observed also in the context of spin systems [ISIIBTI. 

As noted before, the limiting product measures z/^ oo ( 1 1 ) exist for all (j) < cq 
and for reasons explained below, we call the family of measures {z^<^^(p),oo : P > 
O} the fluid phase. The pressure of the fluid phase is given by 

Pfiuid{4>) ■■= lim y log 2:^(0) = log2;oo(0) = log ^° (23) 
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Figure 2: Properties of fluid and grand-canonical measures for cq = 2, ci = 1 as given 



in (23 1 to (26l. Left: Pressure Pgcan (full red line), pjiuid (broken red line) and logZR 
for L = 2, 4, 8 (dashed blue lines), demonstrating the fast convergence to Pgcan- Right: 
Entropy densities Sgcan (full red line) and s fiuu (broken red line). 

and we define the entropy density by the negative Legendre transform 

Sfluidip) := - sup (p \0g4> - Pfluid{4>)) =p(0oo(p)) -plog0oo(p) = 

<P>0 

= (1 + p) log(l + p) - p(log Co + log p) , (24) 



where the supremum is attained for (f) = 0oo(p) ( [T3| ). Note that the fluid pressure 
and entropy density are different from the grand-canonical quantities, because 
^r{4>) = oo for (j) > Cl This yields 

p,„„W:^.imJ.og4W^ -^-J (25) 

and the negative Legendre transform of the pressure is given by 

V.n(p) = | . /^r^^^ ^P<Pc _ (26) 

gcanKHJ | Sfi^id{Pc) " (p " Pc) log Ci , p > p^ 

Note that the Legendre transform of the pressure is usually called the free energy 
density. In thermodynamics, the free energy F is related to the entropy S via 
F = U — TS, where U is the internal energy and T the temperature. Since there 
is no energy and temperature in our case, we define the entropy density as the 



negative free energy density. The functions (23) to (26) are illustrated in Figure 
|2j In analogy to previous results [fTTl [T6l we expect a condensation transition 
for p > Pc. But the non-standard behaviour of the grand-canonical measures, in 



particular the lack of a law of large numbers (22), will lead to a richer behaviour 
than in previous studies, which can be fully understood only in the context of the 
equivalence of ensembles. In particular, the grand-canonical approach alone does 
not provide a complete picture of the phase transition. 
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3.2 Canonical measures 

The canonical measures are given by 

^L,N := 4A • I = AT) where ^^(77) ■= Vx , (27) 



L 



i.e. they are given by a grand-canonical measure conditioned on a fixed number 
N of particles. Their mass functions are independent of (f) and given in terms of 
the stationary weights (|6]) by 

^lAv) = ^ ^^(r?) N) , (28) 

concentrating on configurations 

Xl,n = {r]eXL\ = N] . (29) 

The partition function is now given by the finite sum 

In the following we analyze the limiting behaviour of this quantity. In the dis- 
cussion configurations with many particles on a small number of sites turn out to 
play an important role. Therefore we define the disjoint sets of configurations 

Xl',n = ^ Xl^n \ r]x> R for exactly m sites x G A^} (31) 

with more than R particles on exactly m sites. 

Theorem 1 Suppose R ^ log L, i.e. as L ^ 00. Then the limit 

Scan{p) ■= lim \ log ^L,iv , wheve N/L ^ p , (32) 
exists and is called the canonical entropy density. It is given by 

g (^p-j = Sfluidip) , P < Ptrans ^^3) 

™" \ SfluidiPc) + ScondiP, Pc) , P > Ptrans ' 

where 

Scond{p,Pc) = lim ylogWRUp - Pc)L) . (34) 
The transition density ptrans{o) is given by the unique solution of 

a = {sfiuid{Pc) - (p - Pc) logci - Sfiuidip)^ ^ log ^ , (35) 
where ptrans{(i) > Pc + CL with equality if and only if a = 0. 
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Note that with (|34]) and Q the contribution of the condensate to the canonical 
entropy is given by 

Scond{p,Pc) = -(p - Pc) logci - a log — . (36) 

Cl 



As a special case, taking a = we have ptrans = Pc as the unique solution of p5) , 
and comparing ([33]) with ([26]) yields 



Scanip) = Sgcan{p) for all p > . (37) 

On the other hand, both entropies are different whenever a > 0. 

Proof of Theorem^ Using ( [sij ) we decompose the state space Xl^n = Um=o ^T,n- 
The maximal number M of sites containing more than R particles is certainly 
bounded by M =: \N/R\ . Notice that M p/a E (0, cx)] as L — > oo, and in 
particular M/L ^ 0. We can estimate the number of "uncondensed" configura- 
tions where no site has more than R particles by the following Lemma, which is 
proved in the appendix. 

Lemma 1 For all L,N > 1 and M as above we have 

<|X°^|<|X^,^|. (38) 



This includes for all p >0 and N/ L ^ p 

lim y log \XIj^\ = lim y log \XlM = x{p) , (39) 

L— »oo ij L— voo ij 

where xip) ■= (1 + p) log(l + p) - plogp. 



Furthermore, ifR^ yL, then lim |X° jy|/|X2,,Ar| = \ for all p > 0. 
Now we split the partition function accordingly 

M 

Zl,n = Y.^In. where Z^^ = w^^iX^^) . (40) 

m=0 

For the term m = we get with Lemma [T] 

^ \ogZl^ = ^ log (co^|X°,^|) ^ x{p) - P log Co = Sfi^M . (41) 
The contributions of the other terms are given by 

^'^~\ml^ ^0 \^L-ra,N-k\\ I • l^^J 

^ ^ k=m{R+l) ^ ^ 
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Here we have chosen m sites on which we distribute k particles such that each 
site contains at least R + 1 particles, giving rise to the first and last combinatorial 
factor. The N — k remaining particles are distributed on L — m sites such that 
none contains more than R particles. The sum can be approximated by an integral 
and evaluated by the saddle point method. The saddle point equation reads 

Co L .fN — k\ , k — niR—l ^ ^^^^ 

This has a solution if and only if 

N — {L — m)pc> m(R + 1) or, equivalently p > pc + ma . (44) 
In this case, to leading order the solution to ( [43] ) is given by 

k c:^ N - {L -m)pc (45) 
where we have used that m/L ^ for all m < M. On the other hand, for 



p < Pc + ma the sum in (42) is maximized for the boundary value k = m{R+l) 



Z'^^c^l (^)(|)™Co'^e(^-'")^(^-"^'^) ,ma<p<p,+ma. (46) 



We get in leading exponential order 

, p < ma 

'L,N — \ \mJ \ ciJ ""0 

Q (|)-«+(^-™)p=cr^e(^-™)^(''^) , p>p, + ma 
For p > pc + Q- we get a rough estimate by adding both cases, 

where C = exp ^(co/ci)'''=e~^*^'^'=^ j . Now, to leading order 

— ^log — + ^logL — alog— <0 as L — »• oo , (48) 

L ci L ci 

since M/L ^ 0, R/ L ^ a > 0. This holds only if M -C L/ logL or, equiva- 
lently, R ^ log L. Since p > Pc + a the first summand on the right-hand side of 




(47) vanishes with an analogous argument. Therefore 

M 



^log (i + J2zT,n/zu)^o 

^ m=2 ^ 



(49) 
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and the only exponential contribution to (47) is given by Zj^j^. Thus we have, 
using (|46l), 



1 1 
lim -\ogJ2zZN= j^ogZ 

-L— >oo l7 ' ' — >f^ ' • 



m=l 



L— »00 



L— >00 



lim — log ( L ( — ) 



1 



L-l,(L-l)p<; 



Cl 



(a + Pc) log — - p log Cl + x{Pc) = S fluidiPc) + Scondip, Pc) ■ (50) 
Co 



This is a linear function in p with the same slope — logci as Sgcanip) (26). Note 



that for p ^ oo the first term ( [41] ) behaves as 

Sfiuidip) ^ -P log Co + log(l + p) + 1 . 



(51) 



Therefore, whereas for small p (41) dominates the partition function, (50) dom- 
inates for large p, since it has larger asymptotic slope — logci > — logco- The 
transition density ptrans as a function of a is found by equating both contributions 
which leads directly to ( |35] ). Differentiating the right-hand side of this equation 
yields 



a'(p) = 1 - log ^^-^ / log — 

p / Cl 



(52) 



Thus a'(pc) = and a'(p) G (0, 1) for all p > pc. Since also a(pc) = 0, (35) has 
a unique solution ptransia) > Pc for all a > 0. Further we have 



Ptr 



> 1 for all p> Pc 



a'iPtransia)) 

and thus ptrans{a) > Pc + a with equality if and only if a = 0. 



(53) 
□ 



Now, if a > then M as defined after (31 ) is bounded and converges to p/a, and 
thus (|47]) imphes that 



M 



m=2 



m / rv\ 

L,N ^L,N 



as L — oo . 



(54) 



This is significantly stronger than (48) and it is easy to see that it still holds for 
a = 0, as long as i? ^ vTlogX. Thus for L — > oo the canonical measure 
concentrates on certain parts of the state space, and from the proof of Theorem 
[T (47) the rate of convergence is faster than polynomial in L. Therefore we can 
immediately deduce the following. 
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Corollary 1 For R > log L we have 

P < Ptrans =^ ^L,n{XI^^) 1 , U'-Kl,n{Xl,N \ -^L,Ar) ~^ , 

P > Ptrans =^ ^L.A^ (^L,Ar) ^ 1 , L'^T' L,N {X L,N \ Xl,n) ~^ ) (^5) 

for all n E as L ^ oo and N/ L p. 

This implies in analogy to p3] ), that for p > p trans a typical configuration consists 
of a homogeneous background with density pc and the (p — pc)-^ excess particles 
concentrate in a single lattice site. We expect this kind of behaviour actually 
already for R ^ log L, since wr has an exponential tail and maximal fluctuations 
under in the occupation number are of order logL. Our estimates are not 
strong enough to deduce this, but we are primarily interested in a > 0, which is 
covered by the above result. The same is true for the last statement of Lemma [1} 
For p = Ptrans the contributions of condensed and fluid configurations to the 



canonical entropy are equal (33). This is true on the exponential scale and to 
deduce the behaviour on the transition line we need a finer estimate, given in the 
following Theorem. 

Theorem 2 For N/ L ptrans cind a> Qwe have 

y^R{XlN)/wR{XlN) = 0{L^'^) - oo a. L ^ oo , (56) 
which implies t:l,n{,X\^^) — > 1. 

So in case of a discontinuous transition (i.e. a > 0) the transition line belongs to 
the condensed phase C / F. For a = the transition is continuous and therefore 
p = Pc belongs to the fluid phase F{E). 



Proof. According to (42 ) in the proof of Theorem [T] 



N 

^R\^L,N) — ^^G q \^L-l,N~k\ 



k=R+l 



where k = N—{L—l)pc+o{L) is the solution to the saddle point equation (43 1. In 



addition to the proof of Theorem [T] we consider the next order of the expansion to 
get the correct asymptotic behaviour. Since with Lemma 



_ P(l+P) 

for all p > and k E {R + 1, N), the asymptotic behaviour of the Gaussian 



T)x"(p) 
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integral with variance = — L/x"(pc) (1 + o(l)) is given by its normalization 
and we get 

= L^/2(^)''^"^'"^cr^|X°_,,^(,_,)|v/27rp.(l + pJ(l + o(l)) . (58) 

VCq/ 

With C = ^y2Tcpc{l + Pc) this leads to 

^r{^L,n) _r^T3/2 \co) ^1 \^L-1ML-1)\ , ^/.nN _ 

= CL^''^-^ ^'L'-' ^ 1 + 0(1)) , (59) 

•^0 I L-l 

where we have used the third statement of Lemma [T] that holds for a > 0. We 
use Stirling's formula for the binomial coefficients and note that due to Theorem 
[T]the exponential terms in the ratio vanish, which leaves us with 

"^fiyo'"^! =CL^/^(l + o(l)) ^cx) asL^oo,iV/L->pi,,„,. (60) 
Together with Theorem [T] this implies that 

W^{Xl^N \ ^i,Ar) ^0 as L ^ OO, N/L ptrans , (61) 

which implies the last statement of the Theorem. □ 



4 Equivalence of ensembles 
4.1 Specific relative entropy 

In Table [T] we summarize the results of the previous section in connection with 
the phase diagram shown in Figure [T] In particular, for a = the phases F and 
F/C are empty since pc = ptrans, and we have Scanip) = Sgcanip) for all p > as 



noted already in (37). This implies that the canonical entropy density is concave 
and the condensation transition is continuous. On the other hand, for a > 
we have equivalence of ensembles only in phase F(E), the canonical entropy 
density is non-concave, and the transition is discontinuous. These results concern 
equivalence of ensembles in terms of convergence of entropies of the canonical 
and the grand-canonical measure. In Figure |3] they are illustrated by numerical 
calculations of the canonical entropy density using the recursion relation 

N 

Zl,n = y^^WR{k)ZL-i,N-k ■ (62) 

fc=0 
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phase 


canonical entropy 

Scanip) 


grand-canonical entropy 

^gcan (p) 


F(E) 


Sfluid{p) 


Sfluidip) 


F,F/C 
C/F 


Sfluid{p) 

Sfiuid{Pc)-{p-Pc) log ci-alog ^ 


S fluid{pc)-{p-pc)^OgCi 



Table 1: Summary of the results of Section 3: Comparison between canonical and grand- 
canonical entropy density. Equivalence of ensembles holds only in phase F{E). 




Figure 3: Canonical entropy density Scan{p) for various values of cq, ci and a. Data 
points are calculated numerically according to (62 1 with L = 100 (x), 200 (+), 400 (O), 
and show good agreement with the theoretical predictions for the thermodynamic limit 
(see Table [1]). 
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As can be seen, the grand-canonical entropy density is equal to the concave hull 
of Scanip) which itsclf is not concave for a > 0. The canonical entropy density 
further coincides with the one of the fluid phase, up to the point when it becomes 
metastable and the condensed phase becomes stable. This point has been derived 
exactly by studying the dominating terms in the canonical partition function. 

We can make a connection to other formulations of the equivalence of ensem- 
bles, using the specific relative entropy 

MvTL,^, 4,r) ■= T H{^L,N. 4,r) = 7 ( log ■ (63) 

With the identity t^Y,n — ni -I^l = N), this can be expressed in two useful 
forms, 

h{^L,N, 4,r) = -\ logz/J^(S,. = N) = 

N 1 

= log ZR{(t)) - — log - - log Zl,n ■ (64) 



The derivation of these expressions is straightforward, see e.g. UTTfl . The follow- 
ing is a direct consequence of our results on the canonical measure in Theorem 

m 



Corollary 2 Choosing = 4>r{p) according to (16) we get for all p > 

h{7:L,N, ^ia(p),R) ^ Sgcan{p) " Scan(p) • (65) 



Proof. We use the second expression in (64) for the specific relative entropy. 
Choosing = 0i?(p), the first two terms converge 

log Zii{(I)r{p)) - log (j)R{p) — Sgcanip) (66) 

to the grand-canonical entropy density (26), since with Proposition [1] analogous 
to ([17]) and ^ 

Convergence of the third term in ( 64 ) has been shown in Theorem [TJ which fin- 
ishes the proof. □ 

We can read from Table [T] that 

, P< Pc 

Sgcanip) -Scanip) = { S fiuidiPc)- S fluidip)- (p- Pc) log Ci , Pc<P<PtraM>'&) 

alog(co/ci) ,P>Ptrans 
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In particular, for a = we have pc = ptrans and 
Ktil^n. ^ij,{p),R) ^0 for all p > , 



(69) 



whereas for a > this holds only for p < p^. By a standard result [[61, conver- 
gence in specific relative entropy implies weak convergence, i.e. convergence of 
expectations of bounded cylinder test functions / G Co,b(X), 



as L ^ oo , N/ L ^ p . 



(70) 



This is another formulation of the equivalence of ensembles. 

Furthermore, we can compare the canonical measures with the expected fluid 
measures for the background. 



Theorem 3 Let a > 0. Choosing 



)(p) according to {13) we get 



h{7TL,N, i^^^(p),oo) ^ Sfiuidip) - Scanip) = far < p < Ptr 



whereas far = ci 

Now let a = and R ^ V L log L. We have pt 
< p< Pc, ^far p > Pc- 



{p - Pc) log — > far p > ptran 
Cl 



(71) 



(72) 



Pc and { 71 ) holds far 



Proof. According to the definition ( 63 1 we have 



-—log 

L cq L 



. ^ = Pfluid{(Poo{p)) 

j\ogZL,N + \ 5Z ^L,N{'n)\ogw^{r]) , (73) 



where we have used the definitions ( [TT] ) and ( [28| ), 

1 



w^irj)5{J:Liri),N). (74) 



Splitting the last term of ( [73] ) and using Corollary [T] we see that 

7 Yl ^L,Ni'n) log + J ^^ttlAti) log w^{r]) -p log Co ,(75) 



as L — > oo, N/L p, as long as p < ptrans- Therefore, with definition (24), 

h{7TL,N, ^i^{p),oo) Sfluidip) - Scanip) = for p < ptrans ■ (76) 
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This holds for all a > as long as R ^ \JL log L. For p > ptrans we also use 
(73 ) where 0oo(p) is replaced by Ci. Again with Corollary [T] the main contribution 



to the last term comes now from r) G X\ ^. The sum can be computed by the 
saddle point method analogous to the proof of Theorem [T] and we get 

Y ttl^at (77) log w^(77) ^ a log— -pc log Co - (p-pc) log Cl .(77) 

The same also holds for p = ptrans, since with Theorem [2] vr^ at concentrates on 



Xj ^ also in this case. The first terms in (73 1 are now 



Pfluid{ci) - ^ log — SfiuidiPc) + P log Co - (p - Pc) log Cl , (78) 
L Co 

and together with the behaviour of Scan from Theorem [T] we get for p > ptrans 

H^L,N, l^too) ^ (P - Pc) log - > , (79) 

Cl 

finishing the proof of Theorem [3] Note that a = is included as a special case in 
the above derivation as long as i? ^ ^/L\ogL. □ 



dTTl) allows us to identify the hmit measure and we have 

(/).,^(,,» asL^oo, N/L^p. (80) 

As a direct consequence of the relative entropy inequality ((5], Lemma 3.1), this 
holds not only for bounded cylinder test functions /, but for the larger class with 
(e''-^)iy^ , , < 00 for some e > 0. Since the fluid measures have finite expo- 
nential moments, this includes local occupation numbers f{r]) = rj^, which are 
unbounded. This ensures convergence of densities for p < ptrans ip < Pc for 
a = 0), i.e. in the fluid phases F{E), F and F/C. 

4.2 The condensed phase 



(72 ) may suggest that the limiting distribution of the background in the condensed 
phase C/F is more complicated than the expected fluid measure ^^,00- Together 
with Corollary [T] we can show that the non-zero specific relative entropy is only 
due to the contribution of the single condensate site and indeed the background 
distribution is as expected. In the following we attach some (arbitrary) ordering 
to the lattice sites and identify = {1, . . . , L}. On Xl^i we define the measure 
vTj^ AT as a marginal on the first L — 1 coordinates 

vrL,7V := t!-l,n{-\L G argmax)^'"'-^"^ (81) 
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where 7Tl,n{-\L G argmax) denotes the measure ttl^n conditioned on the event 
that r]L > r]x for all X = 1, . . . , L — 1. Since tcl^^ is invariant under site permuta- 
tions, we have 

vtl.tv := t^l,n{- \y e argmax)^^\^^> for all y E Al ■ (82) 

Note that tt^ tv concentrates on a subset of X^^i, 

Xl-1 ■■= {r] e Xl-1 I Sl-i(^) < N, fi^, <N- ^L-iifj)} , (83) 

and in case of condensation it can be interpreted as the distribution of the back- 
ground. 

Theorem 4 For p > ptrans we have as L ^ oo, N / L p, 

H{7rL^N,i^t^)^0, (84) 
and thus for bounded cylinder test functions 

- ■ (85) 



Note that the first statement ( 84 1 involves the total rather than the specific relative 
entropy and is therefore much stronger than Corollary [2] and Theorem|3j This im- 
plies convergence in total variation norm [5|. Such a result is not possible below 
criticality, since the conditioning on the particle number in the canonical mea- 
sures leads to divergence of the relative entropy. Above criticality, this condition 
is accounted for purely by the condensate site and does not affect the background, 
which shows the same fluctuations as i.i.d. random variables. Following recent re- 
sults in [24], this enables to show that the stationary density profiles converge to a 
Brownian motion with a jump at the location of the condensate. Below criticality 
the corresponding expected behaviour would be a Brownian bridge, but there is 
no proof so far. 



The second statement ( 85 1 is a direct consequence of the first but not a very strong 
one, since it would also follow from convergence in specific relative entropy. The 
site with maximum occupation number will be in the support of the cylinder test 
function only with probability of order 1/L. But due to this possibility, the test 
function has to be bounded, not necessarily by a constant but by a number of or- 
der o(L). This excludes f{rj) = 7]^ as expected, since the expected density does 



not converge for p > ptrans- Note also that with f72] ) and ( [85] ) this system is an 
example where weak convergence is strictly weaker than convergence in specific 
relative entropy. 



Proof. ( 8 1 1 and ( 83 ) imply that 



t<'l,n{^ G argmax) ^^-^ 
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where {fj,N — Si_i(r))) G Xl,n denotes the concatenated configuration. By 
permutation invariance we get 

71l,n{L G argmax) = -^vrL,Ar(X^ jv) + Rl,n = ^ + ^(1)) ^^'^^ 

where < Rl,n < t^l,n{Xl,n \ X\^n)- So the error is exponentially small in the 
system size for p > ptrans (see Corollary [Tj) and of order L"^/^ for p = ptrans- 
Now we can compute the relative entropy 



'^l,n{L G argmax) v^^^^[r]) 



^t^lAvi) log 



(Cl/Co)^^-l(^)Zi,jv (1 + 0(1)) 



+Rl,n , (88) 



where analogous to (31 1 

xl_, = {fieXL-i\vi,-..,VL-i<R} ■ 

Therefore we have 



(89) 



R 



L,N\ 



7^L,N{f])\og 



V, 



L-1 



r)) (1 + 0(1)) 



< C7ri,jv (Xz.,^ \ (X°^^ U X\ ,,)) as L ^ oo 



< 



(90) 



using Corollary [T] and Theorem [2| since the argument of the logarithm is at most 
exponential in L. On X^_i we have w^^^{f]) = c. 



-SL_l(f,) 





and thus 



' Cl,OOj 



-lAXU log (^;/^o)^cr^^I. 4-^ ^ ^ (91^ 



^L,iV [Co — Cij 



where we have used :n-i_Ar(X°_j^) = nL,N{Xlj^). With Theorems 

have for p > ptrans 

ZL,N = w^^{XlA{l + oil)) = 



and 2 we 



Vcq/ 



L-l,pc{L-l) 



1)1 v/27rp,(l + p,)(l + o(l)) ,(92) 



and according to Lemma [T] 

'(l + p,)(L-l)-l 
L-2 



I I 
|^L-l,Pc(L-l)l 



(1 + 0(1)). 



(93) 



A careful application of Stirling's formula, which was not necessary in the proof 
of Theorem |2} yields 



:i+p,)(L-i)-i 

L-2 



Co 



Cq — Ci 



L-1 



c^M (27rp,(l+pe)L)"'/' (1+0(1)) , (94) 
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where we have used in the exponential temi that pc = ci/(co — ci). Plugging 
everything into ( [9T| ) this leads to a perfect cancellation and we get as L ^ oo 

.L-l ■ 



mL,N, = TIlAXIn) log (1 + 0(1)) + 0(1) 



(95) 



which finishes the proof of the first statement. 

Let / G Co,b(X) be a cylinder test function bounded by C and supported on 
the lattice sites supp{f) C N with \supp{f)\ = n. In the following let L > 
max supp{f) such that supp{f) C Al- We have 



where due to Corollary [T] 



(96) 



\R\ 



L,N\ 



^^.^(^) /(^) ^ C7rLAXL,N \ ^ , (97) 



Since |argmax(77) | = 1 for all rj E X\ ^, we have 

Y ^L,N{'n) firj) = J^Yl Yl ^L,N{'n I argmax = {y}) f{r]) 

tYY vrL,^(^ I argmax ={y})^^\^^^/(r7) + <^ 
L — n 



_ 1 
~ L 



Jj,N 



(98) 



due to (82), where boundedness of / implies 



I r2 I 

\-^l,n\ 



J Y Y ^i,Jv(^ I argmax = {y}) /(r)) 



yesupp(f) rjeXl 



< j^Cttl,n{XIA ^ asL-^oo. 



(99) 



Note that this is the only place where we crucially require that / is bounded by a 
constant of order o(L). With Corollary [T] we get 

Y ^LAv)fiv) = {fK. + Rl 



?3 



(100) 



where 



I R3 I 

\-'^l,n\ 



Y ^LAv)f{v) <C7,lAXl^i\xIA 

C nL,N {Xl,n \ {XI n U XlA) ^0 as L ^ oo . 



(101) 
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Together with ([84]) shown above, this implies 

'\f)-'^L,N = {f)^L^N + Rl,n + Rl,n + R\,n {f)uc^,^ > (102) 
since convergence in total relative entropy implies weak convergence. □ 

5 Metastability 

In the previous section the role of the density Pc + a remains open. A first hint 
appears in the proof of Theorem[T] where the saddle point equation for condensate 
contributions ( [43| ) only has solutions for p > Pc + a. But in the context of the 
equivalence of ensembles we are not able to distinguish the phases F and F/ C 
in the phase diagram (Figure [T]), as can be seen in Table [T] A further analysis 
of the canonical measures in terms of the order parameter of the model, i.e. the 
background density pi,g of uncondensed particles, will clarify this point. Since for 
a = the condensation transition is continuous, we only consider the case a > 
throughout this section. We define the observable 

S^^(r7) :=SL(r7)-maxr/, , (103) 

which can be interpreted as the number of particles in the background, since at 
most one site contributes to the condensate. 



Theorem 5 Let Si, S2, ■ ■ ■ E N be any sequence with Sl/L ^ pbg > 0. Then the 
limit 

IM := - lim y \og7iL,N{^L' = Sl) E [0, 00] (104) 

exists for all p > (N/L p), and defines the rate function for the events 
[T,^f = For pbg > p, Ip{pbg) = 00 and for p^g < p it can be written as 

Ipipbg) = Scanip) " S fluid{Pbg) + 

^ f (P - pbg) log Co , pbg> p-a 

\{p- Pbg) log ci + a log(co/ci) , Pbg< p-a 

Proof. For p^g > p, Sl > N eventually and thus 'kl,n {^'l = ^l) = eventually. 
For Pbg < p we use the identity 

ln^R(■^{^L = N}) 
-..» = W«(-|S. = iV)= = iV) 



and the fact that (|64]) and (|65]) imply 
1 



logi'0fl(p),i?(Si=^) = h{7CL,N, ^in{p),R) ^ Sgcanip) " S^anip) ■ (107) 
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Furthermore, Corollary [T] implies that 

lim y log7iL,N{-) = lim 7 log7ri,^( . n {Xl^ U Xl^)) , (108) 
and therefore we get 



1 

L— >oo L 

1 



lim - \0g7lL,N{^L = = Sgcanip) " Scan{p) + 



+ lim - logz/0^(p),R(r7L = N - Sl) + 
+ jim ^ log z/^-^) = r^i, .., r/i_i < (iV - A R) . (109) 

For the last two terms we have fixed the maximum to be on site L, since the cor- 
responding polynomial correction vanishes on the logarithmic scale in the limit. 
With the definition of the single site measure ([8]) the second last term is given by 

lim y log i'4,^{p),R{r]L = N - Sl) = 



(P - Pbg) log(0gcan(p)/Co) , Pbg > P - a 

(P - Pbg) log(0gcan(p)/Cl) - alog(Co/Ci) , Pbg < P - a 



(110) 



where (cf. (fT6l)) 



0oo(p) = Cop/(l + p) , p< p. 

l'-^'^ 1 ci , P> p. 



bgUp) ■■= Jim Mp) = { "^^"^ ' ' ■ (111) 



Due to the condition r/i, . . . , r^^^i < (A^ — S*/,) A i? in the last term, which follows 
from ( |108| ) and ( |103| ), all configurations in that event have the same probability 
and we get 

^im ^ log u^^-l^^^{T.L-i = Sl, Vi, < {N - Sl) A R) = 

L^^L yzRiMp))"--' ' ' 

= Pbg log ^Scanip) _ p(0^^^^(p)) + ^(p^^) = 
Co 

= {Pbg - P) log 0<;can(p) " Sgcan{p) + S fluid{Pbg) ■ (112) 

where 

^Li,s, = e Xl-i,s, hi, ... , r/L-1 <{N- Sl) A i?} . (113) 

Due to the more restrictive condition this is only a subset of Xl_-^ g^, but com- 
pletely analogously to Lemma [T] one can show that as L ^ oo 

^ log \^l~i,Sl \ (l+Pbj?) log(l+Pfc9) - PbalogPfes = Xipbg) ■ (114) 



L 



23 




Figure 4: The rate function Ip{phg) for cq = 2, ci = 1, a = 0.5 and various values of p. 
For p > pc + a the function has a local minimum at phg = Pc, which becomes the global 
minimum for p > ptrans- 



Inserting (110) and (112) into ( 109 ) finishes the proof. 



□ 



Figure [4] shows that the distribution of concentrates on values of the order pL 
for p < Ptrans ciud ou valucs of the order p^L for p > ptrans- These two cases 
correspond to the phases F/ C and C / F, respectively, and have been identified 
already in the previous section. But in Figure |4] also the role oi pc + a can be 
identified. For p < pc + a, the rate function Ip{pbg) ( 104) has only one minimum 
Ip(p) = 0, whereas for p > Pc + a it has an additional local minimum at p^g = Pc, 
i.e. the condensed phase becomes metastable. For p > ptrans this local minimum 
becomes the global one, and the fluid phase becomes metastable. For p = ptrans 
the rate function vanishes for both phases, but the finer analysis of Theorem |2] 
reveals that the fluid phase is already metastable in this case. 

By definition, the observable ^^^(ry) changes at most by ±1 during each jump 
of a particle. So the process (p'f^{r]{t)))^^^ is a one-dimensional simple random 
walk (or a birth-death process) on {0,1,..., A^}, whose stationary large deviation 
rate function is I p. The minima of this rate function correspond to the fluid phase 
for Phg = p and the condensed phase for pbg = Pc- For finite L the system has two 
quasi- stationary distributions 



VrL,Af(-|^LAr) and TTL,Ni-\Xl 



L,N) 



(115) 



corresponding to the fluid and the condensed phase, respectively. Analogous to 



(81 ) we define 



'.\X\j^, L E argmax)"^' 



,L-1 



(116) 



Proposition 3 In the limit L ^ 00, N/L p, we have for all p > 

^(7rL,7v(-|^L,7v).'^0oo(p),oo) ^ ^L,n{-\XI^) l^^^(p),oo , (HV) 
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and for all p > pc + a 



L-1 



' ci,oo 



) — and TrL,Ni-\Xl 



V, 



c\ ,00 



(118) 



In both cases the second convergence is weakly with respect to bounded cylinder 
test functions. 



As in Theorem |4} we can show convergence in total relative entropy ( 118| ) for the 
condensed phase, which is much stronger than convergence in specific relative 
entropy (see comments in the previous section). 



Proof. The first statements in ( 1 17 1 and (118) can be proved analogous to Theo- 
rem|3]and Theorem]?} respectively. Since 



N) 



„ .Ivo irf) 



7; rlyO iri) 

\^l,n\ 



is the uniform measure on X° ^, we get for (117) analogous to (73 ) 

/i(7rL,7v(.|X]^,Ar),Z^^^(p),oo) = 



(119) 



7 7rL,jv(?7|X^,7v)log 



)(0oo(p))' 



Ap)/coT\X[ 
1 



L.N 



I 

l,n\ 



loK X 



P/«md(0oo(p)j - -r log Y 

Li Cq Li 

^ s fluidip) + p log Co - x(/o) = as L ^ 00 



I 

l,n\ 



(120) 



for all p > 0, using Lemma [T[ Note that here we are a priori restricted to X° ^ so 
that there is no error term as in (TTSl). 



The same holds for a modification of (88 ) to derive (118). Using 



TiL Ar(argmax = L, 



largmax=L,Xi ^ (^7, N-lli-lij))) 



L,N) 



argmax=L,X 



1 (r7,iV-Ei_i(T7))(121) 



and 7ri jv(-^i,Ar) — > 1, we get in direct analogy to the proof of Theorem |4] 

L(co/ci)-^cr^zoo(ci)^-i 



H(ttl,n, ^, 



' Cl,OOj 



7rL,N{f])\og 



N 



^L,7v(Xi^^)log(l + o(l)) -^0 asL^oo. (122) 



The second statements in ( |1 17| ) and ( |118 ) follow completely analogously to 
the proofs of Theorem |3] and Theorem |4j □ 
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With Theorems |3] and |4] both statements of Proposition [3] follow directly from 
Corollary [T] and Theorem |2| but only for p < ptrans and p > ptrans, respectively, 
where the quasi-stationary distributions converge to the stationary distribution. 

For finite L, both phases have life-times of the order ~ e^^'^^^ exponential in 
the system size for all p > pc + a, where the exponential rate ^(p) depends on the 
density. It can be calculated using the hitting times 

rf ™'^(p) := inf {t > I maxr/^(t) > R} , 

rl^nd^p-^ ._ |t > I maxr]^{t) < R} , (123) 

which depend on the initial configuration as well as the time evolution. Due to 
the effective one-dimensional random walk picture mentioned above, the quasi- 
stationary expectations of these random variables are determined by the rate func- 
tions at the locations of local minima and maxima. These are 

Ip{p) = Scanip) - Sfiuidip) (min.) 
Ip{p -a) = Scanip) - Sfiuidip - a) + a log Co (max.) 

IpiPc) = Scanip) - SfiuidiPc) + o log — + (p - pc) log Ci (min.) , (124) 

where the last two are only defined for p > pc + a (cf. Figure |4]). Note that 

Ipip) = for p < Pc, whereas Ipipc) = for p > pc. For p > pc + a we then 
have 

e-'^'ip) := jim i log « -^'^(p) »^^_^(.|^o^) /,(p -a)- /,(p) = 

= Sfiuidip) - Sfiuidip - (^) + a log Co , 
r"'(p) := I log « rr'ip) ».,,,(.|xi,,),e^*= Ipip -a)- /p(pe) = 
= Sfiuidipc) - Sfiuidip -a) + ipc + a- p) logci (125) 

where <C . . n{-\x1 f^),e^^ denotes the average with respect to a quasi-stationary 
initial distribution and the time evolution given by the generator C Q. Note that 
for p < Pc + a, Sfiuidip) = oo and ^condip) is not defined, since the condensed 
phase is not stable. The asymptotic behaviour as p — oo is given by 

^flmd^p^ ^ log ^ 

1 + p — a 

^cond/ N ^ ^iQg _ iQgQ ^p_^\^ const. oo . (126) 

Cl 

For all a > we have ^■^'™'^(pc + a) > = ^™"'^(pc + a) and ^■^'™'^(pt 
C,'^"^^iptrans), ^s cxpcctcd. TMs bchaviour is illustrated in Figure [sjfor SOliic spe- 
cific values of the parameters. The predictions are in very good agreement with 



:rans ) 

some spe- 
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Figure 5: Life-times of fluid and condensed phase for cq = 2, ci = 1, a = 0.5. 
Top: The exponential rate ^(p) of the life-time as a function of the density, x and o 
denote Monte Carlo data, errors are of the size of the symbols. Bottom left: Expected 
life-times as used in ( 125 I in a logarithmic plot as a function of L for p = ptrans- Bottom 
right: Tail distribution of the normalized lifetimes r£""'^(pirans)/ T'j^^'^{ptrans) 
compared with the tail of an Exp{l) random variable. 
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data from Monte Carlo simulations, a few of which are presented in the figure. On 
the bottom left for p = ptrans we see that the expected lifetimes for the condensed 
phase are larger than for the fluid phase, which is in accordance with Theorem [2} 
There appears to be a polynomial correction in the condensed phase, but the data 
are not good enough to measure the power in L. 

Note that the last part of the derivation in this section is not rigorous, since 
strictly speaking {tI'I {r]{t))) is not a Markov process. Still one could use a 
potential theoretic approach analogous to [131, to show rigorously that the aver- 
age life times of both phases are exponential in L. However, getting the right 
timescale with this approach would require quite some technical effort. Besides 
the exponential growth rate of the life times with the system size L, simulations 
also indicate that the distribution of the lifetimes is actually exponential, as can 
be seen in Figure |5] on the bottom right. This is to be expected, since the system 
effectively jumps between the two metastable phases in a Markovian way. 



6 Dependence on the number of particles 

In this section we consider the case where the jump rates depend on the number of 
particles in the system rather than the lattice size, which is also the case in some 
models for granular clustering [|23l l4ll. one of our main motivations for this study. 
We modify our original model ([2]), 

^^(^) = {ci ! A^>^ forfc>l, ^?(0)=0, (127) 

where R is now a function of the number of particles SL(r7). For simplicity we 
concentrate on the specific choice R = aJ^L^t}) with a G [0, 1), since a > 1 is 
not interesting for this model. So in principle, the jump rates do not only depend 
on the local occupation number but on the global configuration. But restricted to 
a subset tv with fixed particle number Sl(t7) = N, Ris just a parameter, the 
process is well defined and standard results on stationary measures apply. There- 



fore the canonical measures are well defined as in (28). In particular. Theorems 



[T] to |4] still hold and the proofs apply directly, where a should be replaced by 
ap, since now — > pa. So analogous to (35) the transition density ptrans is 
determined by the relation 

a= (^SfiuidiPc) - (p - Pc)'^ogci - Sfiuidipij ^ (^P^og'^^ , (128) 



and s fluid is given as in (24). The canonical entropy density Scan is still given by 
([33]), but the contribution of the condensate, which determines the behaviour for 
large p, is now given by 

Scond{p,pc) = -p(a\og'^ + logcij +pc log Ci . (129) 
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This leads to 

(•^ fluid (P) ; P — Ptrans 

S fluidiPc) -p(a^Og^+ log Ci^+PclogCi , p> Ptrans ' ^^^^^ 

To study the equivalence of ensembles, one has to define the grand-canonical 
measures. This is not as straightforward as in ([8]), since the number of particles 
Si (77) and thus R is now a random variable. However, we know that the set of all 
stationary measures is convex, and the extremal points are the canonical measures 
(see e.g. [|22ll or |fT6l ). So the grand-canonical measures can be defined as convex 
combinations of canonical measures. 



If the weights WR{k) depended only on the system size L, this would be equiv- 
alent to ([8]), but here the measures are obviously not of product form since the 
weights depend on the total number of particles through R = aT^L^rj). Also the 
normalizing partition function 

00 

Zr{(P):=Y,<P''Zl,n (132) 

Af=0 

does not factorize, since now Z^^^ = w^^{Xl^^). Nevertheless we can define 
the pressure 

^ 00 

Pgcan{(t)) ■= Jim y log ^ (t^^Zi^N , (133) 



N=0 



and by a saddle point argument analogous to the proof of Theorem [T] this is well 
defined and given by 

Pgcani(p) = SUp {p log + Scan{p)) , (134) 
p>0 



the Legendre transform of the negative canonical entropy density (130). With 
( |129| ) we have for p > ptrans 

p\og(t) + Scan{p) = p(^log</)- alog^ - logCij + p^logCi , (135) 

which, analogously to ([25]), implies 
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Co=2, Ci = l, a-0.2 

.«gciin(P) 




0.5 Pmeli,Pc(a) 2 Ptrans 3 
P 



Figure 6: Pressure and entropies for a = 0.2, cq = 2 and ci = 1 as given in ( 136 1, ( 130 1 
and ( |139 l. Data points are calculated numerically according to (62i with L = 100 (x), 
200 (+), 400 {O), and show good agreement with the theoretical predictions for the 
thermodynamic limit. 



The difference is that now the pressure is finite up to 

'Coy 



(137) 



which is strictly bigger than the value ci in (25) for all a G (0, 1). Note that for 



density defined as in ( 14 1 is now a-dependent and given by 

Poo(0c(a)) 



(f) = 0c(a) the saddle point argument (134) does not apply and (133) diverges, 
so Pgcan{<Pcio.)) = oo, scc Figure |6 left. As a consequence of ( 136), the critical 



Pc[a) 



— C 



1-a ' 



(138) 



where poo is still given by ( 12 ). So analogous to ( 26 ), the grand-canonical entropy 



density is given by the negative Legendre transform of ( 136 ), 



^gcan 



''fluid 



Sfiuidip) , P < Pc{a) 

(Pc(a)) - (P - Pc(a)) log 0c(a) , P > Pc(a) 



(139) 



By definition, this is again the concave hull of Scan(p), as can be seen in Figure 
|6| right. The canonical entropy density is calculated numerically using (62) for 
different values of L and A^, and as before the results agree very well with the 
predictions. 

As in the original model, the case a = leads to a continuous phase transition 
and this line of the phase diagram is identical to Figure [T] But for all a E (0, 1), 
Pc(a) > Pc(0), which is the value in ( 14) for the original model. So the phase re- 
gion F(E) in the phase diagram is larger than in the original model (see Figure|7]). 
To complete the phase diagram, we have to derive the analogue of the transition 
line Pc + a, which we call Pmeta in the following. This is defined by the emergence 
of a metastable condensed phase, characterized by a second local maximum of the 
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Figure 7: Stationary phase diagram of the process ( 127 ) for cq = 2, ci = 1. The phases 
F{E), F/C and C/F are defined in Section 2, F(.Eyand F/C overlap (shaded region) 
and phase F is empty. 

rate function Ip{pbg) in Theorem [5} The proof of this theorem makes use of the 
grand-canonical measures, and since these are now of different form, it does not 
apply directly. However, with some effort the proof can be written purely in terms 



of canonical measures (not shown here), and so the result (105) still applies, of 
course with a replaced by p a. An analysis similar to Section 5 reveals that the 
rate function has an additional local minimum 

Ip{pM) forp>p„eta = ^. (140) 

So there exists a metastable condensed phase with background density Pc(0), 
which is still the same as in the previous model, independent of a. This is to 
be expected, since the outflow of the condensate site has to match the background 
current. A simple heuristic argument along these lines provides a general frame- 
work to understand the transition, and is presented in detail in [il8il . Note that in 
comparison with ( 138| ), 



praeta{a) < pc{a) for aU a E [0, 1) , (141) 

with equality if and only if a = 0. This follows immediately from the elementary 
inequality — 1 < (1 — a) (x — 1). So the phase regions F(E) and F/C as 
defined in Section 2 overlap (see shaded region in Figure |7]), and the region F is 
empty. In contrast to our previous model, the equivalence of ensembles still holds 
in the presence of a metastable condensed phase. 



7 Discussion 

7.1 Differences to previous results 

In the following we discuss differences in the condensation transition between 
zero-range processes with and without size-dependence in the jump rates. To 



31 



simplify matters we concentrate on the rates ([2]) for L-dependent jump rates, but 
the features we discuss should hold in general. 

• Without L-dependence the condensation transition in zero-range processes 
is continuous, i.e. the background density pi,g is a continuous function of the 
total particle density p. For model ([2]) this is only true if a = 0, for a > the 
background density p^g = pc < Ptrans is smaller than the transition density 
and the transition is discontinuous. 

• If the jump rates do not depend on L, the equivalence of ensembles holds 
for all densities, and for p > Pc the entropy density is linear in p which is 
often characterized as partial equivalence of ensembles fSHSTl. The reason 
is that the contribution of the condensate to the entropy density vanishes as 
L ^ oo. In model (|2]) this contribution does not vanish, cf. Theorem[T] and 
therefore we have only equivalence of ensembles for p < Pc and nonequiv- 
alence for larger densities. As a consequence of this, the canonical entropy 
density is non-concave, whereas it is concave in case of no L-dependence. 

• Another striking feature of model (|2]) is that it exhibits ergodicity breaking, 
i.e. iox p > Pc + a there are two phases, fluid and condensed, with life- 
times exponential in L, one of which is metastable depending on the density. 
Without L-dependence in the jump rates this does not occur, and for all 
densities there is only one stable phase, either fluid for p < pc or condensed 
for p > Pc. 

So far a discontinuous transition in a zero-range process has only been observed 
heuristically in a two-species system where the stationary state is not known [fTSl . 
The above features only concern the stationary measure, and for systems without 
L-dependence they have been shown rigorously in a general context |fT6l . In the 
following we comment on further differences regarding equilibration and station- 
ary dynamics, which have been studied only heuristically so far. 

• If we prepare a system without L-dependence in a homogeneous distribu- 
tion with density p > it exhibits coarsening [ITTl [T3l . Initially, clusters 
form all over the lattice, and as time progresses the larger cluster sites gain 
particles on the expense of the smaller cluster, leading to a self-similar time 
evolution. The driving force for this behaviour is the fact that there is no 
stable fluid phase with density p > pc. This is not the case in model Q, 
which does not exhibit coarsening for that reason. Instead, it takes a time 
of order e^^'"''*^ before the condensate appears. 

• In a similar setting metastability has been reported as a precursor of the 
coarsening regime, i.e. before coarsening to a single condensate sets in 
[|20ll . Unlike in the present case, heuristic theoretical analysis supported by 
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Monte-Carlo simulation shows that the life time of these metastable config- 
urations does not grow exponentially with system size. 



• For systems without L-dependence, in the condensed phase the distribution 
of the homogeneous background has a sub-exponential tail [|T6ll . In con- 
nection to this, the stationary time scale for movement of the condensate 
location (once a single condensate has build up) is also sub-exponential in 
L, as was found heuristically in [14] in case of a power law. For model 
the background distribution is just z/^^ oo(^) ~ (ci/co)^*^ (see Theorem 
which has an exponential tail. Therefore condensates can move only by 
dissolving completely and, after a time of order e^^'""*^ in the coexisting 
fluid phase, forming on a different site. So the time scale for the stationary 
motion of a condensate is exponential in the system size. 



The long time it takes to form a condensate in the present model is observed in 
Monte Carlo simulations and is explained heuristically by a random walk picture 
in Section 5. The time it takes for the transition between the phases depends on 
the specific model as well as the definition of the phases. In any case its order is 
subexponential in the system size, and for the model ([2]) it is actually of order L. 
Moreover, if p > pc + na for n > 1 also more than one condensate is possible. 
But as can be seen in the proof of Theorem [T] the contribution to the partition 
function of such a configuration is negligible. Therefore one typically observes 
only one condensate, which is a common feature with the stationary behaviour 
of a system without L-dependence, although both cases have very different life 
times. In [|27ll a hydrodynamic theory is developed for the time evolution under 
Eulerian scaling above the condensation threshold. This leads to a generic picture 
for the evolution of a space-dependent initial density profile with total supercriti- 
cal density in systems with rates that do not depend on L. It would be interesting 
to study this problem in the present model. 

Finally, we would also like to stress an intriguing difference to the usual theory 
of first order phase transitions in statistical mechanics. In systems with finite local 
state space or with bounded Hamiltonians, such as spin systems (Ising model) 
or exclusion models, the pressure p is defined for all fugacities > 0, and a 
first order phase transition is a result of the pressure being non-analytic (see e.g. 
[|25l[35ll ). In the model we studied here, the pressure (25 1 is defined only for cp < 
ci, but is analytic on its domain. Therefore the phase transition is a result of this 
bounded domain in connection with the conservation of the particle number, and 
cannot be understood by studying the grand-canonical measures alone. This is in 
contrast to previously studied systems without L-dependent jump rates, where the 
presence of condensation can be characterized by a closed domain of the pressure 
P (cf. nil). 
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7.2 Summary 



In this paper we presented a rigorous analysis of a discontinuous phase transition 
in a simple zero-range process with size-dependent jump rates. The model acts 
as a prototype for systems with that feature and the results are expected to be 
qualitatively similar for a large class of models. This will be discussed in detail 
in a forthcoming publication Going beyond earlier heuristic discussions 

of the phase transition in terms of the order parameter f23[ IH [71, our analysis 
provides a detailed picture of the phase transition in terms of the entire ensemble. 
In particular we note that our approach is a pure equilibrium description, based as 
in the work of ||23l l4| on the notion of thermally activated processes. Hence there 
is no need for an appeal f?| to a non-equilibrium dissipative structure, maintained 
by a flux of entropy, for understanding the nature of the condensation transition 
in the granular shaking experiment. 

Since our results only concern the stationary distribution they do not depend 
on the geometry or dimension of the lattice. The model shows the same features 
observed in granular clustering, namely metastability and a first order transition. 
The only difference is that there is no region in the phase diagram where the fluid 
phase becomes unstable. This is due to the simple choice of rates in this first 
analysis, and the issue will be addressed in [fTSl . 
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Appendix 

A. 1 Proof of Proposition |2] 



For p < pc(22) follows by standard results and for p > p^we. have as L ^ oo 



.l^,,M>H)^^^^{ir. (142) 
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Therefore if we define the truncated occupation numbers fj^ = rjx ^ -R' we have 



= 1 - (l - ^i.(p),K(^. > ^))'^ < CL[^y . (143) 

With R ^ log L this bound is summable and the Borel-Cantelli Lemma implies 
that 

\Y.'^--\Y.^-^^ as L oo. (144) 

Moreover (r},) = + 0((|)«/2) and Var{f,^,) < + 0(i?(|)^/2) and 

therefore by the usual strong law we have J2x&Al ~^ Pc ci-^- ■ 



Taken together, this implies (22) for p > pc, and p = Pc works analogously with 



the power R/2 replaced by _R/4. 
A.2 Proof of Lemma [1] 

Each configuration in jv ^ has at least one site with more than R particles 
and we denote the number of such sites by 

Eiri) := ^ lr,^>Riri) • (145) 

Note that for rj G X^^n \ n ^^^^ 

1 < E{r]) <M= \N/R] , (146) 



where M is as defined in (31 ). For each configuration we define 

S{r}) := [t]^ a R I xe^L) U [ri^-R \ xEAl, Vx>R) G ^L+i?(^),iv • (147) 

If £'(5(77)) > 0, we have to repeat this mapping at most M times such that 

n := S^'iri) e , (148) 

where l{rj) < M denotes the total number of extra coordinates. For Z(r/) < M we 
can identify r) by a configuration in ^, by setting all remaining coordinates 

equal to zero. By this construction it is clear that for each rj G Xl^n \ there 
exists a unique fj E -^l+m.at' l-^- 

\Xl,N \-^L,Ar| = l-^L.A^I " \X^n\ < \Xl+M,n\ ■ (149) 
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Further, each i) has the special property that at least l{ri) sites contain exactly R 
particles, and there are only L sites whose occupation number can be less or equal 
than that. Therefore we can improve the above estimate as 

.0 i^fL + M\., ,fL + M\ \Xl^\ 



\XlM - \xIn\ < [ M )\^U-r\<[ ^ j (L-M)«'^^^^^ 

where the combinatorial factor counts the number of positions of sites with R 
particles. We also used the fact that for all = 1, . . . , i? 

\Xl^^,^,\>{L-M)\Xl^_,\, (151) 

since there are at least L — M positions to put an additional particle without 
violating the constraint 1]^ < R for all x. Together with the obvious fact that 
|-^L,Ar| < |-^L,Af|, this provcs the first statement of the lemma, i.e. 

Xl,n\ < K^L.AT ^ K^L,Af • (152) 



With Stirling's formula we get 



lim ylog(^\^^] = lim f(l+f )log(l+f )-f logf ) =0.(153) 



since M/L^OasL^oo. Therefore 



-log vniiJ^i ^0 (154) 



and (152) certainly includes the second statement of the lemma. More detailed, 
we get to leading order as L — > oo, N/ L ^ p. 

This vanishes for all p > if 

i?logL - ^logi? > logi? V ^ asL^oo, (156) 
R R 

which is certainly the case for R ^ V^- C 



36 



References 

[1] Andjel, E.: Invariant measures for the zero range process. Ann. Probab. 
10(3), 525-547 (1982) 

[2] Angel, A.G., Evans, M.R., Levine, E., Mukamel, D.: Criticality and Con- 
densation in a Non-Conserving Zero Range Process. J. Stat. Mech. P08017 
(2007) 

[3] Bovier, A.: Metastability: a potential theoretic approach. In: Proceedings 
of the ICM 2006, Madrid, pp. 499-518, European Mathematical Society 
(2006) 

[4] Coppex, E, Droz, M., Lipowski, A.: Dynamics of the breakdown of granular 
clusters. Phys. Rev. E 66, 011305 (2002) 

[5] Csiszar, L: I-divergence geometry of probability distributions and minimiza- 
tion problems. Ann. Probab. 3(1), 146-158 (1975) 

[6] Csiszar, I.: Sanov property, generalized i-projection and a conditional limit 
theorem. Ann. Probab. 12, 768-793 (1984) 

[7] Eggers, J.: Sand as Maxwell's Demon. Phys. Rev. Lett. 83(25), 5322-5325 
(1999) 

[8] Ellis, R.S., Haven, K., Turkington, B.: Large deviation principles and com- 
plete equivalence and nonequivalence results for pure and mixed ensembles. 
J. Stat. Phys. 101, 999-1064 (2000) 

[9] Evans, M.R.: Phase transitions in one-dimensional nonequilibrium systems. 
Braz. J. Phys. 30(1), 42-57 (2000) 

[10] Evans, M.R., Hanney, T.: Nonequilibrium statistical mechanics of the zero- 
range process and related models. J. Phys. A: Math. Gen. 38, R195-R239 
(2005) 

[11] Evans, M.R., Hanney, T., Majumdar, S.N.: Interaction driven real-space 
condensation. Phys. Rev. Lett. 97, 010602 (2006) 

[12] Ferrari, PA., Landim, C, Sisko, V.V.: Condensation for a fixed number of 
independent random variables. J. Stat. Phys. 128, 1153-1158 (2007) 

[13] Godreche, C: Dynamics of condensation in zero-range processes. J. Phys. 
A: Math. Gen. 36, 6313-6328 (2003) 

[14] Godreche, C, Luck, J.M.: Dynamics of the condensate in zero-range pro- 
cesses. J. Phys. A: Math. Gen. 38, 7215-7237 (2005) 



37 



[15] Godreche, C: Nonequilibrium phase transition in a non integrable zero- 
range process. J. Phys. A: Math. Gen. 39, 9055-9069 (2006) 

[16] Grosskinsky, S.: Equivalence of ensembles for two-component zero- 
range invariant measures, accepted in Stoch. Proc. Appl. (DOI: 
10.1016/j.spa.2007.09.006) 

[17] Grosskinsky, S., Schtitz, G.M., Spohn, H.: Condensation in the zero range 
process: stationary and dynamical properties. J. Stat. Phys. 113(3/4), 389- 
410 (2003) 

[18] Grosskinsky, S., Schiitz, G.M.: in preparation 

[19] Kafri, Y., Levine, E., Mukamel, D., Schiitz, G.M., Torok, J.: Criterion for 
phase separation in one-dimensional driven systems. Phys. Rev. Lett. 89(3), 
035702 (2002) 

[20] Kaupuzs, J., Mahnke, R., Harris, R.J.: Zero-range model of traffic flow. 
Phys. Rev. E 72(5), 056125 (2005) 

[21] Levine, E., Mukamel, D., Schiitz, G.M.: Zero range process with open 
boudaries. J. Stat. Phys. 120(5/6), 759-778 (2002) 

[22] Liggett, T.M.: Interacting Particle Systems. Springer, Berlin (2004) 

[23] Lipowski, A., Droz, M.: Urn model of separation of sand. Phys. Rev. E 65, 
031307 (2002) 

[24] Loulakis, M., Armendariz, I.: Thermodynamic Limit for the Invariant Mea- 
sures in Supercritical Zero Range Processes. arXiv:0801.2511 (2008) 

[25] Ruelle, D.: Statistical mechanics: rigorous results. W.A. Benjamin, New 
York- Amsterdam (1969) 

[26] Schlichting, H.J., Nordmeier, V.: Strukturen im Sand. KoUektives Verhalten 
und Selbstorganisation bei Granulaten. Math, naturw. Unterricht 496, 323- 
332 (1996) (in German) 

[27] Schiitz, G.M., Harris, R.J.: Hydrodynamics of the zero-range process in the 
condensation regime. J. Stat. Phys. 127(2), 419-430 (2007) 

[28] Schwarzkopf, Y, Evans, M.R., Mukamel, D.: Zero-Range Processes with 
Multiple Condensates: Statics and Dynamics. arXiv:0801.4501 (2008) 

[29] Spitzer, E: Interaction of markov processes. Adv. Math. 5, 246-290 (1970) 

[30] Torok, J.: Analytic study of clustering in shaken granular material using 
zero-range processes. Physica A 355, 374-382 (2005) 



38 



[31] Touchette, H., Ellis, R.S., Turkington, B.: An introduction to the thermo- 
dynamic and macrostate levels of nonequivalent ensembles. Physica A 340, 
138-146 (2004) 

[32] van der Meer, D., van der Weele, J.P., Lohse, D.: Sudden collapse of a 
granular cluster. Phys. Rev. Lett. 88, 174302 (2002) 

[33] van der Meer, D., van der Weele, K., Reimann, P., Lohse, D.: Compart- 
mentalized granular gases: flux model results. J. Stat. Mech.: Theor. Exp. 
P07021 (2007) 

[34] van der Weele, J.P., van der Meer, D., Versluis, M., Lohse, D.: Hysteretic 
clustering in granular gas. Europh. Lett. 53(3), 328-334 (2001) 

[35] Varadhan, S.R.S: Large deviations and applications. Ecole d'Ete de Prob- 
abilites de Saint-Flour XV-XVII. Lecture Notes in Math., Springer, Berlin 
(1988) 



39 



